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SEISMIC PROCESSING METHOD TO IMPROVE SPATIAL RESOLUTION 

This invention relates to the field of seismic data processing. More 
particularly, this invention relates to a method of obtaming improved spatial and/or 
temporal resolution of 2-D or 3-D seismic data. 

A seismic survey represents an attempt to map the subsurface of the earth by 
5 sending acoustic or elastic energy down into the ground and recording the "echoes" 
that return fixmi the rock layers below. The source of the downgoing acoustic or 
elastic energy might come, for example, firom explosions or seismic vibrators on land, 
and air guns in -marine environments. During a seismic survey, the energy source is 
moved across the sur&ce of the earth above a geologic stmcture of interest. Each 

10 time the source is actuated, it generates a seismic signal that travels downward 

Arough the earth, is reflected and/or diffracted, and, upon its return, is recorded at a 
great many locations on the surface. Multiple source-actuation /recording 
combinations are then combined to create a near continuous profile of the subsurface 
that can extend for many kilometers. In a two-dimensional (2-D) seismic survey, the 

1 5 source and recording locations are generally laid out along a single straight line, 

whereas in a three-Hlimensional (3-D) survey the source and recording locations are 
generally distributed across the surface in a grid pattern. In simplest terms, a 2- D 
seismic line can be thought of as giving a cross sectional picture (vertical slice) of the 
earth layers with reflector locations roughly midway between the source positions and 
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the leceiver positions. A 3-D survey produces a data "cube" or voliime ftiat is, at least 
conc^tually, a 3-D picture of the subsurface that lies beneath the survey area with 
reflector positions roughly midway between the source positions and the receiver 
positions in the acquisition grid. 
5 A seismic survey is composed of a very large number of individual seisnuc 

recordings or traces. In a typical 2-D survey, there will usually be several tens of 
thousands of traces, whereas in a 3-D survey the number of individual traces may run 
into the multiple millions of traces. General backgroimd information pertaining to 3- 
D data acquisition and processing may be found in Chapter 6, pages 384-427, of 

10 Seismic Data Processing by Ozdogan Yilmaz^ Society of Exploration Geophysicists, 
1987, the disclosure of which is incorporated h^ein by reference. 

Unprocessed seismic data is only of limited use to an explorationist. Seismic 
data as it is acquired in the field is seldom used directiy, but instead it is first returned 
to a processing center where various mathematical algorithms are applied to the 

IS digital seismic data to enhance its signal content and generally make it more amenable 
to interpretation. A key step in a typical seismic processing sequence is seismic 
imgration. 

As is well known to those skilled in the art, the dip, location and character of a 
reflector on an unmigrated seismic section is rarely representative of the tme dip, 
20 subsurface location and character of tiie structural or stratigraphic feature that gave 
rise to that reflector. Except in the case where the subsurface consists of 
homogenous, horizontal layers, the recorded seismic expression of a stmctural or 



WO02073237 [httpj/AvAV^.getlhepatentcom/Login.dog^^ Page 5 of 31 




WO 02/073237 PCTAJS02/07451 

stratigraphic event must be migrated before it can be reliably used to locate 
subsurface features of interest In areas of steep dip, a reflection that is apparently 
located directly below a particular surface point before migration may, after 
migration, actually be found several hundreds of meters away. Additionally, in 
5 complex stmctural areas where faulting, severe asymmetrical folding and sharp 
synclines are present, dif&actions and multiple reflections may interfere with 
reflections flom the primary reflectors to the point where, without migration, the 
resulting seismic section bears little or no resemblance to the actual subsurface 
stmcture. 

10 Broadly speaking, migration improves a seismic section or volume by 

"focusing" the seismic data contained therein, a process tibat is conceptually similar to 
that of "focusing" the image produced by a slide projector in order to obtain the 
sharpest screen imagp. Migration improves tbe seismic image by correcting the 
lateral mispositioning of dipping seismic reflectors; collapsing diffractions caused by 

15 point scattering centers and subsurface fault terminations; resolving crossing 

reflectors (conflicting dips); and improving the v^cal and lateral resolution of the 
seismic data, among many others. A geneml description of the many ways that 
migration improves seismic data may be found in, for example. Chapters 4 and S, and 
Appendix C, pages 240-383, and 507- 518, of Seismic Data Processing by Ozdogan 

20 Yilmaz, Society of Exploration Geophysicists, 1987, the disclosure of which is 
incorporated herein by reference. As a general rule, seismic data that have been 
properly migrated reveal an enhanced or a truer picture of the subsurface than 
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unmigrated seismic data, the ultimate goal of seismic migration being to produce a 
seismic section or volume that accurately represents flie configuration and character 
of the geology of the subsurface. 

United States Patent 6,049,759 to Etgen discloses a method of piestack 
5 migration of 3-D seismic data. As taught therein, there are two broad variants of 
seismic migration: migration after stack (poststack) and migration before stack 
^restack). Post- stack migration is applied, as the name suggests, to seismic traces 
after they have been stacked—a stacked seismic trace being one that is formed by 
combining together two or more traces after Normal-Moveout (NMO) correction to 

10 form a single composite trace that is an approximation of a zero offset trace at that 
location. Prestack migration^ on the other hand, is applied to seismic traces before 
they have been stacked. Other things being equal prestack migration is always 
preferFed-»both theoretically and in practice-^because it has the potential to produce a 
more accurate picture of the subsurface stratigraphy and structure. However, the 

15 computational effort involved in conq)uting a prestack migration is many times that 
required to do a poststack migratioiL For 2-D seismic hues, this additional 
computational effort is generally manageable and, except for the longest lines, 
prestack migration is often applied to seismic data that has been taken in areas where 
the subsurface structure is complicated. Many 3-D data sets, on the other hand, 

20 contain far too many traces to be cost-effectively migrated via conventional prestack 
algorithms. 

An attractive algoritimi to approximate full prestack migration in 3-D and 

4 
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even for 2-D seismic acquisition is the use of Pre-stack Partial Migration, also known 
as Dip Moveout Correction (DMO). DMO moderates the effect of reflection-point 
smear and enables events to be stacked coherently moderating the effects of dip 
and/or azimuth. Thereafter the data can be stacked and post-stack migrated 
S conventionally but with a substantial degree of data compression. An additional 
advantage is that DMO is a relatively small correction that tends to be insensitive to 
errors in the estimated velocity used for NMO corrections. This is a sig^ficant 
improvement over NMO-stack. 

There are various altOTiative equivalent implementations of DMO. Perhaps 

10 the most popular are Hale's method and the summation method of Deregowski and 

Rocca. Hale's Fourier-based method, proposed in his doctoral Ihesis "Dq> Moveout by 
Fourier Transform" submitted to Stanford University Qeophysics Department, May 
1983, is carried out in frequency/wave vector (f,k) domain. Deregowski and Rocca*s 
summation method described in "Geometrical Optics and Wave Theory of Constant 

15 Offset Sections in Layered Media," Geophysical Prospecting 29, 374406 (1 98 1), is 
carried out in time/space (tpc) domain. It involves summation of data along a "DMO 
trajectory." 

The migration and DMO mediods disclosed in prior art and mentioned above 
start with digitally sampled data from a single seismic line or 2-D grid of seismic lines 
20 and obtain an ou^ut at locations that correspond to nominal positions of locations 
midway between the seismic sources and detectors used in the acquisition process. 
All of the time domain migration and DMO methods involve summation of data 
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moved fix>m a lecoided position to an ou^ut position. In principle this is similar to 
what is performed in NMO correction of the data in a line of seismic data. 

United States Patent 5,596^46 to Wiseciq> teaches a method of preserving 
temporal fiequency components in NMO corrected stack data that are normally lost in 
5 conventional NMO processing. Wisecup teaches the NMO correction of digital 

samples of offset seismic data to zero offset, wh^e the NMO-corrected data may not 
fall at sampling time, and combining all such NMO-corrected contributions prior to 
resampling. This is in contrast to conventional NMO methods wherein an 
interpolation is done on the offset trace prior to NMO correction to ensure that the 

10 NMO-corrected sanq>les &11 on a discrete sampling time. The so-called Random 
Sample Interval Imaging (RSI^^ method of Wisecup preserves more high ftequency 
data than conventional NMO-stack methods. Analogous to and in addition to the loss 
of temporal fiequencies in conventional processing, there is also a loss of spatial 
frequencies in conventional migration and DMO operations. The present invention is 

15 a method of preserving these higher temporal and spatial frequencies in partial or fiill 
migration of seismic data using the concepts taught by Wisecup, 

In one aspect, the present invention is a method for preserving temporal 
frequencies in fiill or partial migration of a line of seismic data using, for example, 
DMO. This differs from the method of Wisecup vx that instead of an NMO operation, 

20 data are migrated to discrete sampling locations defined by the source and/or receiver 
positions. These locations are commonly midway between sources and receivers. 
The migration may be a partial migration such as a DMO, or it may be any full time- 
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or depth-domain migration such as a Kirchofif migration. 

In ano&er aspect of &e invention, the present invention is a method of 
preserving t^poral and spatial frequencies in full or partial migration of a line of 
seismic data. This is similar to the preservation of temporal fiequencies mentioned 
5 above, but additionally, the output spatial locations may be in a location not limited to 
being midway between a source and receiver location. In yet another aspect of the 
invention, only higter spatial fiequencies are preserved. 

Finally, in yet another aspect of the invention, temporal and/or spatial 
frequencies are preserved during a full or partial 3-D migration of seismic data 
1 0 acquired using a 2-D sur&ce grid or a pliuality of 2-D seismic lines wherein the 
ou^ut grid for the migration may be different from the input grid. 

The novel features which are believed to be characteristic of the invention, 
both as to organization and methods of operation, together with the objects and 
advantages thereof, will be better understood from the following detailed description 
15 and the drawings wherein the invention is illustrated by way of example for the 

purpose of illustration and description only and are not intended as a definition of the 
limits of the invention: 

Figs, la, lb show raypaths, travel times and signals corresponding to a point 
diffiractor at a plurality of coincident source-receiver locations. 
20 Fig. 2 illustrates flie migration of an input trace to an output location preserving high 
fiequencies. 

Fig. 3 illustrates raypaths fix)m a point dif&actor located at a position that does not 
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conespond to a surface sampling location. 

Fig. 4 shows fhe results of migration of synthetic data with an input san:q)le interval of 
200 ft to an output sample interval of 5 £t 

Fig. 5 shows the acquisition of a line of marine seismic data for offset souroe- 
S receiver geometry and illustrate a common of&et geometry 

Figs. 6A- 6D show the geometry involved in various migration operatioi^. 
Fig. 7 illustrates a geometry commonly used in acquiring a 3-D seismic data set 

The first embodiment of the present invention is a modification of the 
invention disclosed in United States Patent 5,596,546 to Wisecup and ttie contents of 

10 which are fully incorporated herein by reference. Referring now to Figs, la and lb, a 
schematic illustration of seismic data acquisition using coincident sources and 
receivers is shown. The surface of the earth is indicated by 1 and in the subsurface is 
a point dif&actor 5. The surface locations of the coincident shot-receiver 
combinations are denoted by 10a, 10b, 10c, . . . lOn and the corresponding raypaths 

15 from the surface positions are indicated by 7a, 7b, 7c, . . . 7n. In Fig. la, the vertical 
axis is depth and the horizontal axis is the distance relative to the horizontal distance 
from the point difi&actor. For this isolated dif&actor in the subsurface, seismic traces 
that would be recorded at the sur&ce locations are denoted by 18a, 18b, 18c . . . 18n 
in Fig. lb. In Fig. lb, the vertical axis is time and the horizontal axis is the distance 

20 from ih& point diffractor. The diffraction from the point diffractor are indicated by 

19a, 19b, 19c. . .19n and generally have a time-dependence denoted by fhe difi&action 
curve 15. 
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Denoting by x the horizontal distance between the dif&actor and the recordmg 
position, the time-distance relation of the diflftaction curve for a homograous medium 



Olher, more complex relations are possible for inhomogeneous media. The 
implementation of the first embodiment of the present invention can be better 
understood by referring to Fig. 2. The traces again represent a dif&actor such as 5 in 
Fig. la having an acoustic two-way travel time to the zero offset sensor 6 of 1 ,000 
milliseconds. In a typical seismic survey, acoustic signals 18a • . . 18n at each of the 
sensors 10a flirough lOn are digitized. During digitization of the signals, individual 
samples are taken at spaced apart time intervals of the signals generated by each 
sensor 10a through lOk. The samples are generally shown at 26 for the zero-offset 
sensor 6, and at 24 for the otiier sensors. In flie example of Fig. 2 the samples 26, 24 
typically are taken at a predetermined time interval of two or four milliseconds, and 
are taken at substantially coincident times, that is, all the sensors 10a throu^ lOn are 
simultaneously sampled every four milliseconds. 

During processing of the digitized signals from the plurality of sensors 10a 
through lOn, an ima^ of the dif&actor 5 can be generated by combining the digitized 
signals from each of the sensors 10a through lOn to form a combined output signal. 
As would be known to those versed in the arl, every reflector in the subsurface may 



having a velocity v may generally be given by a relation of the form 




(1) 



9 
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be considered to be made up of a multitude of difBractors, and by proper summing 
along dif&action curves such as 20, a coherent image of a reflector may be obtained 
as a composite of the image of the multitude of difBractors. The diffiaction curve 20 in 
fig. 2 is the same as the curve 15 in Fig. la. In order for the combined signal to 
5 properly image the dif&actor 5, signal combination must occur so that amplitude 
correspondent portions of signals firom each of the sensors 10a through lOn, such as 
the peaks of the wavelets 19a through 19ii, are combined. In generating the combined 
signal, however, amplitude correspondent x>oints, which are shown genmlly at 28 on 
the non-zero offset traces may occur at times other than times at which digital samples 

10 are taken, the sample times being shown generally at 24 as previously discussed. In 
the present embodiment of the invention, all the digital samples 24 in each of the non- 
zero offset traces are time-shifted by a time value which exactly satisfies the 
relationship of eqn. (1) corresponding to the nfiigration distance (x) for each of the 
non- zero offset sensors. 

15 The time values for time-shifting each of the non-zero offset traces are 

calculated to the degree of precision to which a computer (not shown) which performs 
the calculations is capable. Ihe time-shifted samples 24 are stored in a buffer, which 
can form part of the computer memoiy. The combination in the buffer of the samples 
26 firom the zero-offset trace 18e, and the time shifted samples 24 from the non-zero 

20 offset traces forms a combined signal with constructive interference of the diffraction 
signals in a manner similar to that taught in Wisecup, The contents of the buffer 
typically will contain a total number of samples equal to the product of the number of 

10 
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sensois multiplied by the number of samples fiom a single sensor. 

The contents of the buffer also can form a stacked single-sensor signal 
sampled at the equivalent of a much smaller time interval flian the predetermined time 
interval at which the samples 26, 24 were originally taken. 

5 The dif&action curve here and the NMO curve given in Wisecup are similar. 

As disclosed in Wisecup, the NMO stack obtained therein preserves higher temporal 
firequmcies flian prior art methods wherein samples ftom non-zero ofiEset traces are 
interpolated prior to being NMO corrected and sununed. The increase in higher 
temporal frequencies in Wisecup results fix>m the fact that the average sampling rate 

10 can be increased by as much as the multiplicity of traces within the NMO gather and 
fiom eliminating the inherent inaccuracy of interpolation. 

For identical reasons^ the image of the diflfiraction point (and any seismic 
image obtained by a process of time or depth migration) using this first embodiment 
of the invention preserves higher temporal frequency than prior art methods wherein 

15 samples for a non-zero migration distance are temporally interpolated prior to being 
moved along the difEraction curve and summed. As disclosed in Wisecup^ even if the 
input data are aliased in time, frequencies above the Nyquist frequency corresponding 
to the input time-digitization interval are preserved in the output. In an optional 
embodiment of the present invention, the imaged data are low-pass filtered with a 

20 fiiequency ofMx fyyq where /vyg is the input Nyquist frequency and M is an integer 
less than or equal to N where N is flie number of traces in the aperture of the 
migration operator. M is related to the output sample interval with the output sample 

11 
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interval being as small as 1/N of the input sample interval. 

In much of seismic processing and interpretation, obtaining hi^er spatial 
resolution is more of a problem than is obtaining higjier temporal resolution. 
Accordingly, in ano&er embodiment of the invention, the migration is performed with 
5 only the improved spatial resolution, i.e., the diffraction operator is defined with 
respect to close spatial locations at discrete sampling values and the contributions 
firom the non-zero of&ets are time-interpolated prior to being time shifted. 
Alternatively, both improved temporal and spatial resolution may be obtained. 

The improved spatial sampling is illustrated for an acquisition geometry of 

10 depicted in Fig. 3. This is similar to Fig. la in that there is a single diffractor 105 and 
plurality of coincident source-receiver locations 110a, 110b • . .1101, llOj . . • llOn 
with corresponding diffraction raypaths 107a, 107b • • • 107n. A significant 
difference is that the surface location 106 corresponding to a zero migration distance 
does not have an acquisition point located thereat This is a more realistic model. 

15 The discussion above with reference to temporal sampling applies equalfy to the 
spatial sample position represented by the point 106. 

Turning now to Fig. 4, the results of migrating a set of synthetic data are 
shown. The data were generated (not shown) for a model with a point diffractor and 
sampled with a digitization interval of 4 ms. and a spatial sampling of 200 ft. The 

20 traces 151a, 151b, 151c... show the results of nugration of the data wi& the 200 ft. 
san^ling intervals at output positions that differ only by 5 ft Loosely applying 
Rayleigh*s criterion (Fundamentals of Optics, Jenkins, F. A., White, H. E., 3rd edition. 



WO02073237 [http:/AvAV^.getmepatentxorn/Login.dog/$sandy.robb1r^ Page 15 of 31 




WO 02/073237 PCT/US02/07451 

1957, p300.)> a wavelet is assumed to be tesolved if tiie anq)litude of the interference 
(partially deconstructed energy out of position) is no greater than 81 .06 percent of the 
tnie wavelet amplitude. Examination of the plots of tiie reconstructed events shows 
that the amplitude appears to satisfy the Rayleigh criterion wifliin 5 to 10 feet of the 
5 correct location, even tihiough the vapnt data was sampled spatially on a 200* grid 
This is verified upon examination of a plot of the peak amplitude (not shown). 

Coixunonly, seismic data are acquired with a plurality of source positions each 
of which is recorded at a plurality receiver positions. In such situations, it is common 
to use the Dip Moveout (DMO) operation, typically carried out on common offset 

10 seismic data as shown in Fig. 5. Shown in tiie figure is a dif&actor 251 and raypaths 
from source positions Si, S2 to the dif&actor 251 and to receiver positions Ri, R2 
corresponding to the source positions, the source-receiver distance being ;c. The 
DMO operation, as described in Deregowski & Roca, takes data fiom a common 
ofket panel and produces a zero-offset seismic section by summing along a DMO 

15 curve that is well defined, just as the migration of a zero-offset section involves 

summation along a diffiaction curve given by eqn. (1). The DMO operation is well 
known and is not discussed hereafter. In one embodiment of the present invention, 
data acquired along a seismic line are sorted into common offset panels, each of 
which is DMO corrected using the DMO operator and the RSI^ method disclosed by 

20 Wisecup. The individual DMO corrected sections may then be stacked. 

Alternatively, a prestack or poststack Kirchoff migration may be performed on 
one or more lines of seismic data; if poststack migration is used, a DMO may be 

13 
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applied to fhe prestack data. United States Patent 4,964, 103 to Johnson discloses a 
mefliod in which ray tracing is perfonned through a 2-D or 3-D model to determine 
travel times fmm any source to any receiver. These travel times may then be used to 
perform a Kirchoff migration tiiat basically combines data along diffiraction curves 
5 determined from the ray-traced travel times. 

Next, application of the method of &e present invention to 3-D migration is 
discussed. Fig. 6a shows a view of a portion of the subsur&ce with a line 303 
indicated thereon. A single dif&actor point 301 is shown vertically below a position 
307 on line 303. For a coincident shot-receiver geometiy, the curve 309 shows a 

10 difi&action curve corresponding to the dif&actor 301 in the direction of the line 303. 

The vertical axis here is shown as time, but it could also be depth. Also shown in Fig. 
6a is a difi&action curve 311 for coincident shot-receiver positions along a direction 
305 perpendicular to tibie line 303. 

Turning now to Fig. 6b, three dififractois 301, 321 and 341 are shown, with 

15 diffractors 321 and 341 being positioned vertically beneath points 327 and 347 on 

lines 323, 343 respectively. As in Pig. 6a, dif&action curves 329, 349 for coincident 
shot receiver gpometiy are shown in tiie inline direction, i.e., along lines 323 and 343 
are shown. These diffraction curves 329 and 349 are similar to the curve 309 for the 
diffractor 301, and, if the seismic velocities in the subsurface do not change, the 

20 curves would be identical. Also shown in Fig. 6b are diffraction curves 330 and 350 
for the difi&actois 321 and 341 for observation points along line 303. 

The diflfractions observed along line 303 are shown again in Fig. 6c as curves 

14 
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309, 330 and 350. Next, Fig* 6d shows the results of migrating the data of Fig. 6c 
using a correct diffraction operator for line 303. The curve 309 is correctly collapsed 
to the image point 301* but the diffinaction curves 330 and 350 are not correctly 
imaged to a single point: instead, a post migration reirmant denoted by the curves 321 ' 
5 and 341' are flie result of this migration. 

The purpose of a 3-D migration of data from a volume of the subsurface is to 
properly collapse out-of-plane diffraction curves such as 330 and 350 to a single 
point This may be accomplished by a full 3-D migration, such as that taught by 
Johnson, 

1 0 Turning next to the more realistic case of non-zero source receiver o£fr ets, one 

embodiment of the present invention is suitable for migration of data acquired in a 
grid of 2-D seismic lines such as 305, 325.. 345. The same method may also be used 
on a 3-D data set. First a DMO is performed along lines such as 305, 325, 345. This 
was discussed above with reference to Fig. 5b. Following the DMO, the data are 

15 stacked along the lines 305, 325, 345 to simulate zero offset geometry and this is then 
followed by a migration of the DMO-stack data. Those versed in the art would 
recognize that the same procedure could also be used for data acquired along lines 
such as 303, 32, 342. 

An example of a commonly used acquisition geometry for 3-D seismic data is 

20 shown in Fig. 7 wherein sources 422 may be placed along seismic lines such as 420a, 
420b . . . 420n wifii receivers 424 deployed along lines 426a, 426b. . .426ii. For such 
a geometry, fiiU 3-D migration may be used. Alternatively, a method such as that 

15 
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described in United States Patent 4»878^04 to Black et al, the contents of which are 
fully incorporated herein by reference, may be used to perform a 3-D DM0 of the 
data. In one embodiment of the present invention, the me&od of Black or other 
suitable method, is used to perform a 3-D DMO witiii enhanced temporal and/or 
spatial resolution using the methods described above. 

In another embodiment of the present invention, a 3-D Kirchoff migration 
with enhanced temporal and/or spatial resolution is performed using any suitable 
Kirchoff migration method and incorporating the RSI^ method herein. 

While the foregoing disclosure is directed to the preferred embodiments of the 
invention, various modifications will be apparent to fliose skilled in the art. It is 
intended tibat all variations within the scope and spirit of the appended claims be 
embraced by the foregoing disclosure. 
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1 1 . A mefhod of preserving spatial frequency components in spatially and 

2 temporally sampled signals exceeding a maximum spatial frequency related to 

3 a predetermined spatial sample interval, said signals corresponding to a 

4 plurality of sensors positioned at spaced apart locations, said method 

5 comprising: 

6 (a) for each of said plurality of sensors, generating time adjusted samples 

7 at an output location, said adjusted samples generated to provide 

8 amplitude correspondence between each of said signals generated by 

9 each of said plurality of sensors using a migration operator related to 

10 different propagation times of energy to each of said plurality of 

1 1 sensors from associated energy sources; and 

12 (b) combining all of said temporally adjusted samples into a combined 

13 output signal having a spatial sampling interval smaller than said 

14 predetermined spatial sample interval. 

1 2. The method of claim 1 wherein said migration operator is selected from a 

2 group consisting of (i) a post-stack time migration operator, (ii) a post stack depth 

3 migration operator, (iii) a pre-stack time migration operator, (iv) a pre-stack depth 

4 migration operator, and, (v) a DMO operator. 

17 
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1 3. The mefliod of claim 1 wherein said sensors and said associated energy 

2 sources define one of (i) a single seismic line, (ii) a grid of seismic lines, and, (iii) a 3- 

3 D seismic data set 

14. A method of preserving spatial and temporal firequency components in 

2 spatially and temporally sampled signals exceeding a maximimi spatial fi:equency and 

3 a maximum temporal frequency related to a predetermined spatial sampling interval 

4 and a temporal sampling interval, said signals corresponding to a plurality of sensors 

5 positioned at spaced apart locations, said method comprising: 

6 (a) for each of said plurality of sensors, generating adjusted samples at an 

7 output location, said time adjusted samples generated to provide 

8 amplitude correspondence between each of said signals generated by 

9 each of said plurality of sensors using a migration operator related to 

10 different propagation times of energy to each of said plurality of 

1 1 sensors firom associated energy sources; and 

12 (b) combining all of said time adjusted samples into a combined 

13 output signal having a spatial sampling time interval smaller than said 

14 predetermined spatial sample interval and a temporal sampling interval 

1 5 smaller than said predetermined temporal sampling interval. 



18 



WO02073237 [httpr/M^vw.getmepatent.conn/Loginxog/SsandyTO Page 21 of 31 




WO 02/073237 PCT/US02/07451 

1 5. The mettiod of claim 4 herein said migration operator is selected from a 

2 group consisting of (i) a post-stack time migration operator, (ii) a post stack depth 

3 migration operator, (iii) a pre-stack time migration operator, (iv) a pre-stack depth 

4 migration operator, and, (v) a DMO operator. 

1 6. The method of claim 4 wherein said sensors and said associated energy 

2 sources define one of (i) a single seismic line, (ii) a grid of seismic lines, and, (iii) a 3 

3 D seismic data set 



1 7- A method ofpreserving at least one frequency component in sampled signals 

2 exceeding a maximum fiiequency related to at least one predetermined sample 

3 interval, said signals corresponding to a plurality of sensors positioned at spaced apart 

4 locations, said method comprising: 

5 (a) for each of said plurality of sensors, generating adjusted samples at an 

6 output location, said adjusted samples generated to provide amplitude 

7 correspondence between each of said signals generated by each of said 

8 plurality of sensors using a migration operator related to different 

9 propagation times of energy to each of said plurality of sensors from 
10 associated energy sources; and 
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11 (b) combining all of said adjusted samples into a combined output signal 

12 having a sampling interval smaller than said predetemiined sample 

13 interval. 

1 8. The method of claim 7 wherein said at least one frequency is selected from the 

2 group consisting of (i) a temporal frequency, and, (ii) a spatial frequency. 

1 9. Hie method of claim 7 wherein said at least one sample interval is selected 

2 from the group consisting of (i) a spatial sample interval, and, (ii) a time sample 

3 interval. 

1 10. The method of claim 7 wherein said migration operator is selected from a 

2 group consisting of (i) a post-stack time migration operator, (ii) a post stack depth 

3 migration operator, (iii) a pre-stack time migration operator, (iv) a pre-stack depth 

4 migration operator, and, (v) a DMO operator. 

1 11. The method of claim 7 wherein said sensors and said associated energy 

2 sources define one of (i) a single seismic line, (ii) a grid of seismic lines, and, (iii) a 3- 

3 D seismic data set. 
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